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EXECUTIVE SUMMARY 


This Technical Note reports the results of reduction 

and preliminary analysis of airborne-laser sea-ice 
surface profiles obtained in April 1977 from Fram 

Strait between northern Greenland and Spitzbergen. 

Data reduction consists of the removal of discontinuities 
in the profiles, and numerical computer filtering to 
remove aircraft motion. Analysis consists of computer 
calculations of descriptive statistics. The data set 
chosen is in proximity to an ice camp (Ruby) operated 
by personnel from the Polar Research Laboratory, Santa 
Barbara, California. Aerial photography taken 
simultaneously with the laser profiles serves to place 
the profile data precisely with respect to Camp Ruby on 
12 April 1977. 
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I. BACKGROUND 


The Polar Oceanography Branch of the Naval Ocean Research and Development 
Activity (NORDA) has on file analog tapes of laser sea-ice surface profiles 
obtained in the Arctic from P-3A aircraft at approximately six month intervals 
since 1970. The total amount of data in the library is about 100,000 nautical 
miles, of which only a small fraction, perhaps 1%, has been reduced and analyzed 
for specific projects. Sea-ice surface profiles provide statistical information 
such as spatial frequency of pressure ridges, distributions of ridge heights and 
spacings, and general sea ice roughness. These parameters are used in either 
numerical models of sea ice movements and distribution, acoustic models in which 
sea ice serves as a boundary or a source of ambient noise, or both. 


The raw data appear on FM analog tapes. Reduction of the data is made 
difficult by the presence of aircraft motion (long-wavelength variations in the 
profile) and various automatic and accidental discontinuities introduced by the 
recording system. Removal of these extraneous effects has been the subject of 
humerous papers (Holyer et al., 1977, and Hibler, 1972). Computer programs have 
been developed which remove aircraft motion numerically, but discontinuities are 
usually of such an unpredictable nature that they must be treated on an individual 
basis by a trained analyst, and cannot be removed automatically by computer. It 
is this individual treatment which makes the data reduction process very slow, and 
has resulted in the large percentage of data which remain unreduced. The Polar 
Oceanography Branch is presently developing the capability to reduce this data 
much more efficiently by means of real-time graphic analysis using a table-top 
minicomputer. This note, however, describes the scheme presently in use. 


TI. DATA REDUCTION SCHEME 


Data produced on an x-y analog plotter directly from an analog laser profile 
tape are shown in Figure 1. The slope in the data between points C and B is one 
exampie of the long-wavelength aircraft motion superimposed on the terrain profile. 
The "event" at point A is a ridge about 4.6 m high. The vertical exaggeration of 
about 30:1 best displays the discontinuities. "Events" at B and C are discontinui- 
ties which must be somehow removed so that the "real" profile can be recovered 
from the data. If these discontinuities are ignored in the data reduction, the 
program which removes aircraft motion will produce spurious peaks at these 
locations. These peaks will be counted as ridges. Point B is an automatic 
shift of the data (back to center-scale) which prevents recordings of analog 
voltages larger than about 1.5 v (a limitation of the recording system). Point C 
is a type of discontinuity called a "dropout". Dropoutr are almost always a total 
change of about one-fourth of full scale, but have widely varying shapes. The 
phenomenon of dropouts is yet to be explained (by anyone, including the manufac- 
turer of the laser profiling system, Spectra Physics). They appear to occur 
predominately over rough terrain. 


Data are usually chosen for reduction on the basis of its proximity to a 
specific experiment area, such as an acoustic array. Location of the data in time 
and space is accomplished by comparing the time (simultaneously recorded in flight 
on a parallel channel of the analog tape) with an aircraft log giving aircraft 
position. Precise determinations of location can then be made by comparing the 
laser profile with aerial photography taken simultaneously with the profile. 


Data are first plotted on an analog recorder and visually reviewed for quality. 
Figure 2 shows typical data hand-annotated with the results of a visual review, 


with one area marked as unsuitable for further reduction because the number of 
discontinuities is considered excessive. "Hand" removal of so many discontinui- 
ties might lead to unreliable descriptive statistics. The "unsuitable" section 
(Fig. 2) where several closely spaced data shifts are shown may indicate that the 
aircraft was executing a maneuver such as an altitude change or turn. This 
condition would result in the laser beam pointing obliquely at the terrain, thus 
giving erroneous ridge heights. 


When it has been decided which profiled area is to be reduced, the profile is 
digitized using an A/D converter, onto a 7-track computer-compatible tape. Data 
are recorded in millivolts (-1500 to +1500 mv) and are sampled at about 200 Hz. 
Since the aircraft ground speed is usually 180 to 250 knots, this amounts to 
an actual data point spacing of 0.46-0.64 m. Average data point spacing used is 
0.5 m. One digital record is usually 15,000 points, or about 7500 m in length. 


Figure 2 also shows time and data~point-number scales added by hand to the 
analog trace. This procedure approximately locates discontinuities in the digital 
record. The approximate locations can then be used to obtain a print listing of 
the small portions of the record near the discontinuity, or the entire record can 
be computer-listed (in millivolts) and the locations used as a guide in finding 
the discontinuities and removing them. 


Table 1 shcws a computer listing of data in the vicinity of a data shift. A 


notation is made to add +1600 mv to each data point after point number 5127. If t 
another discontinuity is reached later in the record, further corrections are then | 
necessary. 


Table 2 shows a computer listing of data in the vicinity of a “dropout.” A 
series of corrections, amounting to subtracting a sloping straight line, is shown 
in Figure 3 (see Holyer et al., 1977). The plot (Fig. 3) of the data shows the 
extent of the discontinuity. 


The cumulative corrections are applied to the data in the computer, and checked 
by tisting and inspecting the corrected record. When all corrections are satis- 
factory, the data are ready for numerical filtering and computation of statistics. 


These two final steps are performed in one computer run by an algorithm sugges“ed 
by Hibler (1972). 


The procedure outlined above has been used on a recently acquired data set to 
demonstrate the output of the computer programs and the presentation of results. ; 
III. DATA NEAR SEA ICE CAMP RUBY, 1977 

In April 1977, personnel from NORDA’s Polar Oceanograp*y Branch operated 
sensors aboard the Naval Oceanographic Office's (NAVOCEANO) P-3A BiRDSEYE aircraft 
as part of a field experiment in the Fram Strait (between Greenland and Spitzbergen). 
Data were gathered around two drifting sea ice stations, Ruby and Pear’, which were 


operated by personnel of the Polar Research Laboratory, Santa Barbara, California. 
The locations of the camps are shown in Figure 4. 


Acoustic experiments were performed in the sea in the vicinity of Camp Ruby, 
and airborne laser profiles and 9 x 9 in aerial photographs were made of the sea 
ice terrain on 12 April. The profiles and photography provided the relative 
locations and statistical terrain features of the sea ice near Camp Ruby for 
12 April, and for as many days before and after 1% April, as the ice remained in 
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nearly the same configuration near the camp. (On 12 April at 1517Z, the camp was 
located at 8217N 0652W, according to the LTN-51 inertial navigation system aboard 
the aircraft). Figure 5 is an uncontrolled photomosaic of the area of Camp Ruby, 
using the 9 X 9 inch cartographic quality aerial photographs. 


Figure 6 is an analog trace of the uncorrected profile data. In this parti- 
cular data set, all data were considered reliable. File and record numbers 
referring to the digital data tape are annotated on the trace, and serve to name 
the record. Many discontinuities (a total of 49) are to be found. Each must be 
removed individ.ally by the procedure described in the previous sections. 


After the data are resampled and numerically filtered with the existing 
FORTRAN program, each final record contains exactly 4980 data points, with aircraft 
motion and discontinuities removed. 


IV. SUMMARY OF RESULTS 


Appendix 1 contains all final reduced data (with terrain height in meters) 
obtained from the uncorrected data plotted in Figure 6. 


Appendix 2 is a reproduction of the computer listings summarizing the data in 
each record in Appendix 1. The ridges counted are listed individually by height 
("Z(M)"), position in the record (data point "I"), and distance from the beginning 
of the record in kilometers ("KM"), which provides scaling to the actual terrain. 
Provision is made for 120 ridges in the listing, but in this particular analysis, 
excess ridges are hand annotated on the page. Excess ridges are located by inspec- 
tion of the appropriate data in Appendix 1. Another table in the listing, headed 
"STATISTICS OF COUNTED RIDGES," gives the height ranges (bins) into which ridges 
are placed according to height, the number in each range, and the number per 
kilometer in each range for this record. A third table headed "PROFILE STATISTICS," 
is based on the entire record, not individual ridges, and contains the mean, root- 
mean-square (RMS), etc., of the data. Other information provided in the listing 
includes the aircraft groundspeed, data point spacing, number of data points in the 
record, actual length of the record in kilometers, number of ridges counted, and 
the location of the highest ridge in the record. 


Ridges are counted by the Rayleigh criterion (horrowed from optics), where a 
local maximum terrain height (over 0.5 m in this case) is counted as a ridge if 
adjacent minima in both directions are less than one half the maximum. No 
minimum ridge spacing or maximum terrain slope is considered (see Diachok, 1975). 


The print listing 1s designed to provide a general analysis of the terrain 
included in the record. It is expected that individual applications of the data 
may use only a small portion of this information. 


Figure 7 shows the location of each record with respect to the terrain. This 
step is performed using either a plot (not shown) of the final reduced data, or an 
analog trace such as Figure 6, and comparing ridges seen in the photography with 
those in the profile. The user can employ Figure 7 to map required ice parameters 
as a function of position. Since, in effect, each data point (terrain height) is 
now located, it is possible to re-analyze the data in any desired way and tie the 
results to an experiment conducted in the area. For one example, one may use a 
sliding window to compute a continuous RMS, or ridge spacing as a function of 
position. 
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Figure 4, Map showing location of Sea lee Camps Ruby and Pearl on 12 April 1977 
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Figure 5. Uncontrolled Photomosaic of sea ice near Camp Ruby on 12 April 1977 
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